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2 ' Abstract 

We study the ttN imN reaction around the A^*(1440) mass-shell energy. Considering the total cross 
sections and invariant mass distributions, we discuss the role of iV*(1440) and its decay processes on this 
reaction. The calculation is performed by extending our previous approach [Phys. Rev. C 69, 025206 (2004)], 
in which only the nucleon and A(1232) were considered as intermediate baryon states. The characteristics 
, in the recent data of the vr^p — > n^n'^n reaction measured by Crystal Ball Collaboration (CBC), can be 

understood as a strong interference between the two decay processes: iV*(1440) ttA and A^*(1440) 



' -^(^^)s%vavc- ^l^*^ found that the scalar-isoscalar tttt rescattering effect in the N N* {■rniY^^^^^ vertex, 

which corresponds to the propagation of a meson, seems to be necessary for explaining the several observables 
\^ , of the ttN ttttN reaction: the large asymmetric shape in the tt'^tt" invariant mass distributions of the 



TT p ^ TT^w°n reaction and the w^p 7r+7r+n total cross section. 
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I. INTRODUCTION 



In many pion induced reactions on the nucleon, the single pion production reaction vrA^ ttttN 
has been studied with a particular interest due to its role as a major inelastic process. The 
partial wave analyses of this reaction, together with other reaction channels such as the elastic and 
ttN rjN reactions, have revealed various properties of the nonstrange baryon resonances, A^*'s 

b eenp er formed 
„ „ „ J 'i' 30,3 and 

i9i,y,yy. 

In Ref. [1^, we have discussed the total cross sections of this reaction up to TV = 400 MeV, 
making use of the chiral reduction formula proposed by Yamagishi and Zahed [ij]. We focused 



and A's (e.g. see Refs 
on the basis of the phenomeno 
the chiral perturbation theory 



. A number of theoretical investigations have also 
ogical approaches using the effective Lagrangian 
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on the role of A(1232), i.e. the influence of vrAA and pNA. interactions on the reaction processes. 
Because these interactions are not directly observed through the two-body decay of A(1232) in 
contrast to the vrA^A interaction, their coupling constants are difficult to determine. We found 
that the vr^p — > tt^tt^p reactions are sensitive to the vrAA and pA^A interactions, and could be a 
source of information of their coupling constants. 

Besides being useful for clarifying the properties of A(1232), the ttN — > vrvrA reaction is expected 
to provide us valuable information about the Roper resonance A* (1440) and its decay to the 
vrvriV channel. iV*(1440) decays to the vrvriV channel via the A^*(1440) vrA and iV*(1440)_-* 
A(7r7r)^^g^^g processes, which have the branching ratios of about 25 % and 7.5 %, respectively 3|- 
The importance of the latter process has already been pointed out in several studies P, Q]. 

The importance of A* (1440) and its subsequent decays has been discussed extensively also in 
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IS 



IS 



20|. However, such two-pion decay 



the pp — > ppiT^-K^ and pn ^(vrTr)^^'' reactions 
of A* (1440) is not easy to study in the elastic vrA ttA scattering, and also in the 7 A vrvrA 
reaction because A* (1440) has small electromagnetic transition rate compared to other relevant 
resonances such as A(1232) and A* (1520) Q- 

Recently, the tt~p — > vr'^vr^n reaction has been measured up to p^- = 750 MeV/c (i.e. T^^ ~ 



620 MeV) by Crystal Ball Collaboration (CBC) 
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These high precision data cover the energy 



region far from the threshold, in particular around the A* (1440) mass-shell energy, Tj^ ~ 480 MeV. 
Several new interesting results related to A* (1440) and its decay processes were reported: (i) The 
energy dependence of the total cross section shows a shoulder-like shape just below the A* (1440) 
energy (see Fig. 14 in Ref. 23 1, and it should be compared with the jp n^iT^n total cross section 
displayed in the same figure, for which no such shape appears because of the small radiative 
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coupling of A^*(1440)). (ii) The Tr^n invariant mass distribution shows a peak near the invariant 
mass equal to the A (1232) energy, which would be produced via the process A^*(1440) — > vrA. (iii) 
The vr'^vr'^ invariant mass distributions shows a large asymmetric shape in population of the events, 
i.e. the peak at large value of m^(7r'^7r^) is larger than that at small value of m^(7r'^7r^). 

When the total energy increases, the correlation between the outgoing pions will become visible. 
The scalar-isoscalar correlation of two pions via the A^*(1440) — > ^{1:71)^^^^^^ decay is particularly 
interesting, because such correlation may generate the a meson pole. In view of the sizable con- 
tribution of iV*(1440) ^(7rvr)^=^^^g, the ttN vrvriV reaction is a possible source of information 
about this controversial scalar meson. It is worth noting that several literatures have suggested that 
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this "o"" degree of freedom is important also in understanding the structure of A^*(1440) 

In this paper, we investigate the vrA^ — > vrvrA^ reaction in the energy region up to = 620 
MeV, especially around the A^*(1440) mass-shell energy. Through the comparison with the recent 
CBC data, we particularly discuss the role of A^*(1440) and its decay processes on this reaction. 
Furthermore, we try to discuss the possibility of extracting the information about the a meson 
such as its existence. The calculation is performed by extending and improving the theoretical 
framework of our previous study ^13] in which the contributions of A(1232) have been discussed 
in detail in the energy region up to T-j^ = 400 MeV. 

This paper is organized as follows. In Sec.^ we give a brief summary of our previous study and 
explain new ingredients introduced in this work. The numerical results are presented in Sec. II I II 
and the contributions of A^*(1440) to the vrA^ — > vrvrA^ reaction are discussed. Then our results 
are compared with the recent CBC data. Summary and conclusions are given in Sec. IIVI In the 
Appendices we summarize some details of phenomenological treatment in the calculation. 

II. THEORETICAL TREATMENT OF nN nnN REACTION 
A. Brief summary of our previous study 

The starting point of our previous study 3] is the chiral Ward identity satisfied by the invariant 
amplitude MmrN- Assigning the four-momentum and isospin indices of the external nucleons and 
pions as in Fig. ^ we have 

M-rrnN = (M^ + M A + M SA + My a) 

+ {ki,a ^ -ks,c) + {k2,b ^ k^.c) 

+Maaa^ (1) 




FIG. 1: The ttN — !■ nnN reaction. Pions have the isospin index (a, 6, c) and the four-momentum ki 
{i = 1, 2, 3), and nucleons have the four-momentum pj [j = 1, 2). 

where the symbol ( <-> ) represents a permutation of the momentum and isospin indices of the pion 
in the first four terms, and 

= ^[{h - k2f - m2]<5-^(Ar(p2)|vr^(0)|iV(pi)), (2) 



= -^(^2 - hrd'^'iNip^m^mNip,)), (s) 

2 /. 

MsA = -^k^S^^ J d^xe-'^''^-''^^ 

x{N{p2)\T*{a{x)j%^m\N{p,)), (4) 

MvA = -i^{ki + k2Ykle'''' j <f xe~'^''^-''^> 

X {N{p2 ) I T* ify^ {x)j%M)\N{pi)), (5) 

Maaa = j^k'^kl^k^ J d^xid^x^e-'^^^^^'^^^^ 

xV(P2)|T*01^(xi)jl,(x2)jl,(0))|iV(pi)). (6) 

The pseudoscalar density is the interpolating pion field with the asymptotic form 7r"(x) 

T^ia^outix) + ■ ■ ■ {xq =Foo). The one-pion reduced axial current is defined by j'X^ix) = 

-|- f-nd^T^"'{x) where A'^{x) is the ordinary axial current with the asymptotic form vl^(x) 
— f'Kd^'K\n,ovLt{x) + • • • {xq Too). The vector current and the scalar density are represented by 
jy^ix) and (t(x), respectively. 

The Ward identity (IT|-(IH1) was first derived by Yamagishi and Zahed making use of the chiral 
2a] . Owing to this formula the invariant amplitude is expressed in terms of 



reduction formula 
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Green's functions of well-defined current and density operators. Then the consequences of broken 
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FIG. 2: The diagrams considered in our previous study Q| to evaluate Green's functions in Eqs. Q-®- 
Crossed versions are also considered for all of these diagrams. The double line represents the propagation 
of the nucleon or A(1232). Note that A(1232) does not propagate in Msa because the iV-A transition is 
not brought about by the scalar density operator a. The pseudoscalar density tt and the vector current jy 
are dominated by the pion and p meson pole, respectively. 

chiral symmetry subject to the asymptotic condition d^^A'^^^x) — > f-nrn^Tii^^outi'^) + • • • {xq — > =Fc>o), 
are exactly embodied on the amplitude without relying on any specific model or expansion scheme. 
Thus, by using the chiral reduction formula at the beginning of the discussion, we can consider 
the detail of each reaction mechanism separately from the general framework required by broken 
chiral symmetry. This nature of the chiral reduction formula has a great advantage in tackling on 
the hadronic processes in the resonance region in which the systematic chiral expansion scheme 
becomes difficult to be implemented. 

Green's functions (i.e. the matrix elements of current and density operators) appearing in 
Eqs. 0-0 are not uniquely determined by broken chiral symmetry. Therefore we need to em- 
ploy a model in order to evaluate Green's functions. In Ref. they were calculated by taking 
a phenomenological approach based on the relativistic tree level diagrams as shown in Fig. [21 
We considered only the nucleon and A(1232} as the intermediate baryons, and vr and p as the 



intermediate mesons. Details of our model 



13| | are summarized in Appendix 1X1 



B. Contribution of iV* (1440) 



Now we extend our previous approach by including A^*(1440). As for the contributions of 
A^*(1440) on the reaction, we consider the ttNN*, N N* (tttt)^^^^^^, and vrAA^* interactions. The 



first two interactions have already been considered in many theoretical investigations of the vrA'^ — s- 
ttttN reaction near threshold. While the vrAA^* interaction has often been neglected, this interaction 
will become important in the energy region considered in this paper. Indeed, the importance of 



N* An decay is suggested experimentally 
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ttNN* and ttAN* interactions 



The ttNN* and vrAA^* vertices with the one-pion leg are generally related to the matrix elements 
of axial current j^^ for the N-N* transition as 



(7) 



and for the A-A^* transition as 



AN* 



Ta(3 1- 



UN*{P), 



(8) 



respectively, where q'^ = {p' — p)^, Q'^ = —q^ and t = (p' — p)'^, t'^ is the isospin Pauli matrix and 
is the j ^ i isospin transition (2i + l) x (2j + l) matrix. The isoquadruplet Rarita-Schwinger 
vector-spinor and the isodoublet Dirac spinor are denoted as U^^{p) and Ui{p) {i = N,N*), respec- 
tively. By employing Eqs. HB5|) and ()B6|) as the effective irNN* and vrAA^* interactions, some of 
the form factors in Eqs. Q-Q can be exactly related to the renormalized coupling constants 



JttNN* [t) — —p- 

JTT 



-F 



NN* 



A,2 



(9) 



= ^ F^fit) + {mr,*-m^)F^f{t) + tF^f{t) 

Jit L 



(10) 



At tree level, all the form factors are reduced to constants. Because it is difficult to fix all of 
their value in the present status of experimental data, we eliminate F^^* ^^^d F^^* by using the 
PCAC hypothesis f-KNN*,-KAN*{'m'^) — f-wNN* ,ttAN*{^)- Then, at tree level, we obtain the analogs 
of Goldberger-Treiman relation for the ttNN interaction. 



f (^2 \ _ pNN* 



(11) 




<y + (\UAjA)s 



FIG. 3: The diagrammatical interpretation for the chiral reduction of the N N* {TnT)^'^^^^^ vertex. The decay 
ampHtude is decomposed into two contributions arising from the different origin in the chiral structure: one 
includes (a) the scalar N*-N transition matrix element, and another includes (b) the N*-N transition matrix 
element of the contact, scalar-isoscalar combination of two axial currents. The former is due to the explicit 
breaking of chiral symmetry and thus vanishes in the chiral limit, whereas the latter does not. 



and 



UAN'{ml) = ^ F^f + {mN*-mA)F^,f , (12) 

In ^ 



respectively. 



Furthermore, in this paper we neglect F^^* foi^ simplicity. This would be partly justified by 
the fact that, in the case of vrA^A interaction, the contribution of the Fa,2 term is considerably 
small compared to the Fa^i term 13]. The form factor F^^' (t) does not appear in our calculation 
as long as we consider the Lagrangian (|B6|) for the ttAA^* interaction. 

Using the central values of the A^*(1440) ttN and N*{1A40) vrA decay widths hsted in 
the particle data table of Ref. Q], we obtain /jrA^Af* (w-^) = 0.465 and fnAN*{m'^) = 1.71. As for 
fnAN*iiT^n)j wc take account of the finite width of A(1232) in the same manner as in Ref. |2l| . 
Those values lead to F^^* = 0.63 and F^^* = 1.15, respectively. 



2. NN*{TnT)jf°^^^ interaction 

We next consider the N N* (tttt)^^^^^^ interaction with the two-pion leg. Recently, we have 
discussed in detail its general and phenomenological aspects on the basis of chiral reduction for- 
mula 27]. We make use of our results also in the present work. 

Using the chiral reduction formula, we find that the N N* {7nr)g^^^^^ vertex is generally described 
by the following matrix elements of currents and density operators (see Fig. One is the scalar 
matrix element of the N-N* transition [Fig. 01 (a)], which is factorized as 

{Nip')\amN*{p)) = -^f^UN{p')uN*{p). (13) 
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N 



N 



(a) 



(b) 



N* 



N* 



FIG. 4: The "minimal model" for the N N* {TnT)^^^^^^ vertex. The form factors cr_Rjv(t) and FAA{t) are 
dominated by the processes: (a) the contact term and (b) the tttt rescattering together with the contact 
scalar- isoscalar baryon-pion interaction (the meshed blob represents the tttt rescattering in /= J=0 channel). 



Another is the matrix element of the N-N* transition caused by the contact, scalar- isoscalar 
combination of two axial currents [Fig. |31 (b)], which is expressed as 



dW^- {N{p')\T*{j%^{x)jUmN*{p)) 
i6''''g^u Faa it)uNip)uN* (p) , 



scalar-contact 



(14) 



where k represents the four-momentum of external pion with the isospin index a. Note that -FyiA(i) 



does not contain an y si ngle baryon poles owing to the definition for the A^*(1440) N{'kit 



decay given in Ref. 



l y si : 



|/=0 

5" wave 



\I=0 



M 



We then obtain the general expression for the N N* {'k-k) ^ ^^^^ 



vertex as 



NN*(tt-k)s 



f2 

J TT 



<yRN{t) —^FAA{t) 



un{p')un* (p). 



(15) 



When aRj\f(t) and FAA{t) are constants, this expression reduces to the result obtained from the 
effective chiral Lagrangian to order 

Instead of taking aRj\f{t) and FAA{t) as constants, we proposed in Ref. ^3] the "minimal model" 
which explicitly includes the scalar-isoscalar correlation of two-pions by considering the tttt rescat- 
tering mechanism in I=J=0 channel (Fig. |1J). This model gives the following parameterizations of 
aRNit) and FAA{t), 



(16) 



FAAit) - Faa X { 1 + lGit)ti=^{t) ] , (17) 



where G{t) and t^-°(t) are the pion loop integral and / = tttt rescattering amplitude, 



' I (27r)4 f-ml + ie {I - Pf - ml + ie ^^^^ 



8 



with P'^ = t, and 



,1=0 us _ 6 t-mi/2 



f^l + il/mt-ml/2)Git)' 
respectively |27l |. Based on the dimensional regularization scheme with a renormalization scale 
/i = 1.2 GeV, the loop integral G{t) can be expressed as 

for i > 4m^, 



= ^^-l + ln-^ + aln^-.vra) (20) 



ml , fj + 1 



for t < 0, and 



2 



= -^^-p- 1^-1 +ln— J +o-(7r - 2arctano-) ) (22) 



for < t < 4ml, where a = - {Aml/t)\. 

Using Eqs. (|15() - H17() and performing the phase space integral, we obtain the following expression 
for the iV*(1440) ^ iV(vr7r)^=^^^g decay width 

rNN'{.n)s = aic*i? + [i{c*2? + ic\c\, (23) 

where c* = —o^pf ji^vruD and = Faa/^. The numerical value of the coefficients a, /?, and 7 
is a = 1.199 X 10-3 GeV^ [3 = 14.06 x lO^^ GeV^ and 7 = 7.754 x lO^^ GeV^ respectively. 
We note that these values are different from those of Refs. Q, This difference arises from 
the vrvr rescattering and the relativistic effects. Both c\ and C2 are not fixed by the A^*(1440) 
N {ttttY^^^^^ decay width, and we can take any value on the ellipse (^5]) on the c[-C2 plane. 

C. Cut-off factor 

We shall introduce the cut-off factors, which stem from the finite size of hadrons, for each vertex 
in a phenomenological manner. In Ref. [l3| we treated all hadrons as the point-like particles. They 
are, however, the bound states of quarks and gluons and thus have the finite size. Although in 
principle those effects should be directly derived from QCD, it is very difficult and still an open 
question which deserves the theoretical challenges. Instead, we consider those effects phenomeno- 
logically according to the discussions of Ref. j2^]. We attach the four-dimensional cut-off factors, 
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(a) - / • • ^ 

(3,0aJa)s 

(b) 

J A J A J A 

^ ^^^ ^ 

(c) 

7.1 J A J A 

X ^^^^ 

^ N* A N 

FIG. 5: The diagrams including iV*(1440) taken into account in the present paper. Ah crossed versions of 
these diagrams are also considered. 



to each leg a of the vertex (where we assume Ua = ^ ) ■ Here pa and rn-Q, are the four- momentum 
and mass of the particle corresponding to the leg a, respectively. The cut-off factor fa is normalized 
to one on the mass-shell, and thus exhibits its effect only for the internal lines. There are five cut-off 
parameters: A^v, Aa, and Ajv* for the intermediate baryons, and A^ and Ap for the pion and p 
meson poles appearing in Figs, andlSfd). In this paper we assume A^r = Aa = A^r* = A^ 
and A^ = Ap = Am- 



III. RESULTS AND DISCUSSIONS 

Putting together all diagrams introduced in our previous study (Fig. [^J and this work (Fig. \^ , 
we calculate the vrA^ — > vrvrA^ total cross sections up to T^r = 620 MeV. We also calculate the vr'^Tr'^ 
and TT^n invariant mass distributions for the tt~p — > vr'^vr'^n reactions whose experimental data 
were reported by CBC [2^ . We especially try to discuss the contributions of A^*(1440) around its 
mass-shell energy. 



A. Fixing parameters 

Before showing our results, we mention how to fix the parameters: the cut-off parameters A^ 
and Am, and the coupling constants of the vrAA and pNA interactions (denoted as fnAA and 
fpNA, respectively), and the parameters and which characterize N N* {'K'k)^^'^^^^^ interaction. 

Because the value of f-^AA and fpNA can not be directly determined by the two-body decay 
width, they are usually extracted from the data using meson exchange model, or estimated by using 
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FIG. 6: The cut-off dependence of the total cross sections. The results are shown for the cases of = 
Am = 750 McV (dashed line), A^ Am = 850 MeV (solid hue), and A^ = Am = 950 McV (dashed-dotted 
line). The results without considering the cut-off factors are shown by the dotted line. The thick solid 
line ill (c) corresponds to the case of Ar = 850 MeV and Am — 2Ar — 1700 MeV. The data are from 

Refs. 00QQ0Q0,QQQQQQQQQ 



references therein). For 



3C|. 



theoretical approach such as the quark model and QCD sum rule. Many reports on their value fall 
into the range of 0.4 < /ttAA ^ 0.8 and 3.5 < /pAfA ^ 7.8 (see and re: 
instance, the quark model relations lead to /ttAA = 0.8 and /pArA = 5.5 

In Ref. Ij], it is found that, at least in the low energy region up to = 400 MeV, the 
influences of vrAA and pA^A interactions are negligible for the 7r~p vr^vr^n and 7r~p vr'^vr^n 
channels. The NN*{'KiT)g^2^^^ interaction (i.e. and C2), on the other hand, gives a negligible 



contribution to other three channels, i.e. vr^p — > tt^tt^p and vr+p vr+vr+n |^. This fact allows 
us to consider {fpNA,fnAA) and (0^,02) separately in the parameter fixing. Thus we first try to 
investigate the influences of Ab, Am, fwAA, and fpNA on the total cross sections of vr^p — > tt^tt^p 
and TT^p TT+TT+n channels, and to estimate the value of these parameters. 



1. Cut- ojf dependence 

In Fig. El we show the cut-off dependence of the total cross sections. In this calculation, we 
include all diagrams but Fig. |S{a) which includes the N N* {Tnr)g^^^^^ vertex, and take /ttAA = 0.4 
and fpNA = 7.8 which are in the range explained above^. 

First of all, if we do not take account of the cut-off factor H24|) for each leg of the vertices, i.e. 
if we treat the hadrons as point-like particle, the resulting total cross sections obviously overshoot 
the experimental data above T.^ ~ 300 MeV. 

^ The reason for taking these value will be discussed later. 
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T„(MeV) T,(MeV) 

FIG. 7: The dependence of numerical results on /^^aa f^nd fpNA- (a) the ir^p tt+tt'^p channel with 
(LaaJpNa) = (0.4,7.8) (solid line), {f„AA,fpNA) = (0.6,7.8) (dashed line), and {J^aaJpNa) = (0.8,7.8) 
(dotted line), and (b) the ir^p — > tt^tt^p channel with (/ttAA, fpNA) = (0.4, 7.8) (solid Hne), (/ttAA, JpNa) = 
(0.4,5.7) (dashed line), and [J-kAAtIpNa) — (0.4,3.5) (dotted line). The contribution of Fig. (Sja) is not 
taken into account. We use Kb — 850 MeV and Am — 1700 MeV. The data are the same as in Fig. IHl 

We choose three typical values = Km = 750, 850, 950 MeV, and calculate the cross 
sections for each value. As shown in Fig. |HJ we can see clear difference among these results above 

~ 300 MeV. For the n^p vr^vr^n channel, the case of A^ = Am = 850 MeV seems most 
appropriate. The difference is obviously seen in the results for the ir'^p tt'^tt^p channel. The 
case of Kb = Km = 950 MeV is quite different from the data compared to the other two cases. For 
the iT~p -K~7r^p channel, all results somewhat underestimate the data. 

Next we consider several cases of Kb^Km- Aj\/ = 0.5Kb, 1.5Kb, 2Kb for each value of A^. As 
for the TT'^p TT'^iT'^n and TT~^p vr+yr'^p channels the results change only within 10 % when we 
vary the value of Km, whereas the Tr~p 7r~-K^p channel is sensitive to the variation of Km- We 
find that the larger value of Am for each A^ seems appropriate: the thick solid line in Fig. IHl^c) for 
Kb = 850 MeV and A^ = 2Kb = 1700 MeV (this line should be compared to the thin solid line 
with Ab = Am = 850 MeV). 

In the following discussions, we take A^ = 850 MeV and Am = 1700 MeV. 

2. The ttAA and pNA interactions 

Here we try to examine how the total cross section changes according to the variation of fnAA 
and fpNA in their allowed range. For the n'^p ■n'^iT'^n channel, the change is small, and the 
results still compatible with the data. For the ir'^p — > tt'^-k^p channel, the influence of the pNA 
interaction can be hardly seen for 3.5 < fpNA ^ 7.8, whereas the ttAA interaction generates visible 
change on this channel [Fig. Efa)]. The lower value of /ttAA seems appropriate to the data. In 
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-10 10 

Ci [(GeV)-'] 



FIG. 8: The ellipse (23) representing the allowed values of and Cj. Here Tj^N*(TrTT)s = 26.25 MeV is 
used, which is the central value listed in the particle data table Q^. The point A, which is {cl,C2) = 
(—4.68,0) GeV^^ in our case, corresponds to the point often used in the study of hadron reactions. The 
points B and C correspond to (c^, c^) = (0, -1.37) GeV^^ and (cj, c^) = (-14.2, 3.92) GeV"\ respectively. 
See the text for the meanings of the solid and dotted parts of this ellipse. 

contrast, for the variation of fpNA the change of the ir^p — > tc^tt^p channel is clearly seen [Fig.ITfb)] 
but is negligible for that of fnAA- The higher value of fpNA seems appropriate to the data. 

This is why we used /ttAA = 0.4 and fpNA = 7.8 in the discussions of cut-off dependence of the 
total cross sections. Also in the following we continue to use these values for /ttAA and fpNA- 



3. The dependence on c* and c\ 

As mentioned in Subsec. IIIBl the values of c{ and C2 can not be fixed by the relation H23|l 
with the decay width of A^*(1440) -^(■^^)s'^avc' which just represents the ellipse on the c\-c\ 
plane as shown in Fig. |HJ Although the case of C2 = and cj < (the point A in Fi g. El ) is 
conventionally used in the phenomenological models of hadron reactions (e.g. see Refs. [sl. l2l|). it 
was first pointed out in Ref. Q] that the contribution should be also considered in view of the 



general chiral effective 
the values of c* and C2 



/agrangian. Several attempts have been recently performed to determine 



30- 

Calculating the total cross sections for all five channels up to = 620 MeV with allowed values 
of c\ and C2, we find that the appropriate values would be on the solid curve in Fig. |H1 The values 
on the dotted-curve generate obvious disagreement in the total cross section, which is consistent 
with other studies of the vr vrvr A^ reaction P] . 

In Fig. |nj we show the total cross sections calculated with several values of and on the 
solid curve in Fig. |H1 The recent CBC data for the 7r~p ir^TT^n channel are plotted as the filled 
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FIG. 9: The c* and C2 dependence of the total cross sections. The sohd hne is the result with (0^,02) at 
the point A on the ellipse, whereas the dashed and dotted lines are at the point B and C, respectively. 
The dotted-dashed line is for (c*, c*) = (—8.0, 1-08) GeV~^. The data are the same as in Fi g. IHl for (a)-(c), 
whereas are taken from Refs. 

00S0 for (d) and (e). The CBC data are plotted as the filled square in (e). 



square in Fig. ini(e). We notice the CBC data exhibit a bump below the iV*(1440) energy, i.e. 350 
<T^< 500 MeV, and this bump is called the "shoulder" in Ref. 

For the ir^p — > tt^tt^p channels, the results show only a small increase in the total cross section 
above T^^ = 500 MeV except for c| and C2 around the point C [see Figs. IHtb) andlHl^c)]. Below 

= 500 MeV, however, the A^A^*(7r7r);g.^g^^ interaction has no influence on the total cross section. 



This result is consistent with other studies 



1^ 



However, the vr p — > tt+vt n and tt n ^ vr'^Tr^n channels are sensitive to the variation of and 
4 [F«. 0d) andETe),. Th. . expected .o. ot.e. ...e. . low ene.g. flS. T.e .e.* 

with and near the points B and C are obviously incompatible with the experimental data. 
The preferable value of d[ seems to be found between the points A and C in view of the data. 

The TT^p TT~^TT^n channel is also sensitive to the variation of and C2 [see Fig. IHta)]. The 
total cross section above T^r = 400 MeV is increased by the contribution of the process Fig. [Sl^a): 
the lower the value of c| becomes, the larger its increase is. Because this channel is, however, 
already reproduced well without A^*(1440), accurate data of the tt~^p — > TT~^TT~^n reaction above 

= 400 MeV will strongly constrain the value of and C2 not to disturb this present status. 
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FIG. 10: The contribution of each diagram including iV*(1440) to the total cross sections. The value of 
(c|,C2) is taken as (—8.0, 1.08) GeV~^. Each line expresses the result including (i) all diagrams (solid line), 
(ii) no TrAiV* vertex (dotted-dashed line), (iii) no N N* {-K-Ky^"^^^^ vertex (2-dotted-dashed line), and (iv) 
neither contribution (dotted line). The data are the same as in Fig.^ 



From the present situation of -7r"'"p ■K'^ir'^n data, we find that c* > —8.0 GeV~^ is preferable in 
our model. 

It is worth noting that, in contrast to other three channels, our results for the ir'p vr^vr^n 
and TT^p — > TT^TT^n channels show peak structure around the A^*(1440) energy corresponding to 

~ 480 MeV, except for those with extremely low values of cj near the point C. This feature 
should be compared with the CBC data for the TT~p ir^-K^n total cross section which shows the 
shoulder at the A^*(1440) energy. 

Considering all above results, we choose = —8.0 GeV~^ (this leads to = 1.08 GeV~^) as a 
plausible value for this parameter. 



B. Each contribution of the diagram including A^*(1440) 

Now we discuss the contribution of each diagram including A^*(1440). The results are shown 
in Fig. 1101 We find that the process Fig. Etb) gives negligible contribution for all channels. The 
process Fig. |Sfc) does not give visible influence on the Tr~^p vr+vr+n and vr+p -K^ir^p channels, 
but brings 30-50% increase to the total cross section above T-j^ = 400 MeV for other three channels. 
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In the iT~p — > TT^TT^n channel, the process including N N* {Tnr)^^2ave vertex dominates the total 
cross sections near the threshold, whereas its contribution decreases gradually above T^r = 400 MeV. 
Instead the contribution of the process Fig. \^c) grows in this energy region. We find that the 



23] can be understood 



shoulder of the n^p -K^-K^n data at the A^*(1440) energy reported in Ref. 
as a result of interference between these two processes rather than the sole contribution of A^*(1440) 
pole. 

With regard to the diagrams including A^*(1440), the vr^p vr+vr^n and ir^p -K^n^n channels 
show similar behavior. It is thus expected that the shoulder at A^*(1440) energy may be observed 
also in the Tr~p — > 'ir~^TT~n channel. To see this, however, more accurate data in those energy region 
are necessary. 

Here we notice that our result underestimates the data above ~ 500 MeV in the ir'p 
ir^ir^p, 7r~p 7r+7r~n, and TT~p — > Tr^ir^n channels (particularly remarkable for n^p ir^n^p). 
This is not surprising because that energy region is above the Roper mass-shell energy. The 
higher resonances such as A^*(1520) and A^*(1535), which do not consider in this work, will become 
relevant to these channels. In particular, A^*(1520) would play a key role since this resonance has 
a large branching ratio about 45 % for the vrvrA^ decay [whereas the A^*(1535) has only about the 
5 % branching ratio]. Indeed, the total cross section data for these three reaction channels show a 
small bump around the A^*(1520) energy (see e.g. Fig. 5(b) and (c) in Ref. Q for the ir'p tt"tt^p 
and TT^p 'K^'K~n channels, and the CBC data 23] for vr^p n^-K^n channel.). However, such 
structure_is not seen in the Tr+p — > 7r+7r"'~n and n'^p vr+vr^p channels (see e.g. Fig. 5(d) and (e) 
in Ref. 



le 1 



C. Invariant mass distribution of tt p ^ n^Tr'^n reaction 

Figure ^2 shows the ir^n invariant mass distributions for several values of T^^. Our result 
captures the qualitative features of the data well. At = 288 MeV, the mass distribution is 
almost determined by the process including N N* {TT7r)g^^^^^ vertex. The contribution of this process 
decreases gradually when becomes higher. Instead, the contribution of the process including 
ttAA^* vertex dominates the mass distributions. As suggested in Ref. |3]) a peak near the A(1232) 
energy, which can be found above T^^ = 488 MeV, is indeed generated by the A^*(1440) — > Avr 
process. Although the A(1232) mass is taken as ttt-a = 1232 MeV, the peak occurs at the energy 
m?{'K^n) ~ (1.210)^ GeV^. In Ref. ^23] the authors said that the shift of the peak is natural 
because the nucleon pole term, which is large background for the ttN elastic channel, has small 
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FIG. 11: The n'^n invariant mass distributions for the ir^p — > Tr'^w'^n reaction for several values of 
(Ptt)- The CBC data are plotted as the histogram. Each line expresses the results including (i) all 
diagrams (solid line), (ii) no ttAN* vertex (dotted-dashed line), (iii) no N N* {7111)^^^^^^ vertex (2-dotted- 
dashed Hne), and (iv) neither contribution (dotted line). The m^(7r°n) — (1.210)^ — 1.464 GeV^ and 
m2(7r"ri) = (1.232)^ 1.518 GeV^ corresponds to the real part of the A(1232) pole and the Brcit-Wigner 
mass of A(1232), respectively. 



contribution to this reaction. However, in view of our numerical result, we find that the peak 
position is obtained as a result of the interference of several processes. 

We next consider the vr'^vr'^ invariant mass distribution (Fig. I12j) . For low T^^, the results are 



dominated by the process including A^A^*(7r7r)^ ^^^^ vertex, which is the same as the case of vr' 



n 



invariant mass distribution. We observe that our result reproduces the large asymmetry in the 
mass distributions around the A^*(1440) energy: a small (large) peak in small (large) value of 
m^(7r'^7r'^) and a depletion in between (see the result of Tj^ = 488 MeV in Fig. I12|) . In our results, 
the process with N N* {TnT)^^^^^^ vertex decreases the distribution for small value of m^(7r*'7r''), 
whereas it increases for large value of m^(7r''7r*^). Without this process, the large peak at large 
value of m^(7r'^7r'^) can not be reproduced. This result indicates that the large asymmetry in 
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FIG. 12: The tt^tt" invariant mass distributions for the tt p ^ Tr'^ir'^n reaction for several values of T^r {Pir)- 
The meaning of each line and histogram is the same as in Fig. 1111 



the tt^tt" mass distribution at the A^*(1440) energy is due to the strong interference between the 
decay processes A^*(1440) Avr and A^*(1440) —>■ N{TnT)g^^^^^. Our numerical results up to 
= 500 MeV agree with the data qualitatively. Above = 500 MeV, however, our results do 
not reproduce a large peak at larger m^(7r'^7r'^). This would be because the contributions of higher 
mass resonances such as A^*(1520) and A^*(1535) are not taken into account. 

It is worth mentioning that the above features of the and TT^n mass distributions at the 
TV* (1 440') energy are characteristic of the two-pion decay of the Roper resonance, A^*(1440) 



Ntttt 



13 



3 



. This indicates that A^*(1440) indeed gives visible effect to the observables of ttN 



ttttN reaction in those energy region. 



D. I = J = TTTT rescattering in NN* {Tnr)g~^^^^ vertex. 

Finally we discuss the tttt rescattering representing the scalar-isoscalar vrvr correlation explicitly 
in the N N* {Tnr)g^^^^^ vertex, which is depicted in Fig. Efb). It has been suggested in several 
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FIG. 13: The tttt rescattering effect in the ttN — > ttttN reaction: the (a) ir^p 7r°7r"n and (b) 
TT+p TT+TT+n total cross sections, and (c) the tt'^tt*' invariant mass distribution of ir^p —>■ Tr'^n'^n 
at — 488 MeV. The dashed hne is the full result including the tttt rescattering with (c| ,02) = 
(—8.0,1.08) GeV~^, whereas the solid line is the result including no tttt rescattering but with the read- 
justed value (cj,c^) = (-13.0, 1.89) GeV-^. 

literatures (see e.g. Refs. 71, [t^) that the a meson pole can be dynamically generated around 
450 — 225i MeV by such rescattering in I=J=0 channel. In view of the vrA^ ttttN reaction 
being sensitive to the process including N N* {Tnr)^^^^^^ vertex, we can expect that this reaction 



also becomes a source of information about this meson. However, no direct evidence of a meson is 

First we try to calculate the tt~p vr'^vr'^n total cross section without the process including tttt 
rescattering in the N N* {Tnr)^^^^^^ vertex, i.e. we consider only the contact interaction Fig. Ufa) 
for this vertex. We find that, even if the tttt rescattering is not considered, we can readjust c| 
and C2 within their allowed values so as to reproduce our results including the vrvr rescattering 
effect^. For instance, the total cross section with = —8.0 [i.e. the solid line in Fig. Efe)] is 
reproduced by choosing = —13.0 and = 1.89 [see Fig. IT^ a)]. The similar results are also 
obtained for the 7r~p vr+vr^n and vr^p tt^tt^p total cross sections, and for the TT^n invariant 
mass distributions of ir^p ir^-K^n. Therefore, although the ^ N* {inr)^^^^^^ vertex is necessary to 
explain the vrA^ ttttN reaction, it seems difficult to conclude whether the a meson pole dominates 
the form factors of the vertices as long as we consider only these total cross sections and invariant 
mass distributions. 

However, the situation is somewhat different in the ir^p vr+vr+n total cross section and the 
vr'^Tr*' invariant mass distributions of vr^p TT^ir^n. The above value (cJ,C2) = (—13.0, 1.89) GeV""*^ 
is not acceptable for the tt~^p — > n'^iT^n total cross section because these values lead to the nu- 



Note that in this case tlie coefficients in the decay width formula H23|l become as a = 0.476 X 10"^ GeV^ 
/? = 5.29 X 10"^ GeV^, and 7 = 2.98 x 10"^ GeV^. Thus the range of allowed values of and Ca changes. 
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merical results far from the data [Fig. EJb)]. Also, it is remarkable that the large asymmetric 
shape in the vr'^vr'^ invariant mass distribution at Tj^ = 488 MeV can not be reproduced without 
considering the vrvr rescattering effect. These results suggest that the vrvr rescattering effect in the 
N N* {ttttY^^^^^ vertex is necessary to our calculation of vrA^ — > ttttN reaction. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the vrA^ ttttN reaction in the energy region up to = 620 MeV, especially 
around the A^*(1440) mass-shell energy. Being motivated by several interesting observations related 
to A^*(1440) in the recent CBC experiment of the ■K~p — > TT^TT^n reaction, we have discussed the role 
of A^*(1440) and its decay processes in the reaction processes. The calculation has been performed 
by extending the theoretical approach constructed in Ref. [l^j ]. 

We have found that A^*(1440) shows a significant contribution to the vrA^ — > vrvrA^ reaction 
through the decay processes iV*(1440) ^ iV(vr7r)^=^^^g and iV*(1440) ^ vrA. In contrast, the N* 
ttN process just gives a negligible contribution. While the contribution of A^*(1440) —>■ N{7nr)g^2^^^ 
process already appears in the threshold region in several channels, the A^*(1440) — > vrA process 
becomes important above = 400 MeV. The characteristics of the CBC data for iT~p -K^Tr^n 
are generated by a strong interference effect between them. 

We have also found that, above = 400 MeV, the tt~^p — > vr+vr+n total cross section is 
remarkably sensitive to the variation of and within their allowed values. Because this reaction 
is almost saturated by the contributions of the nucleon and A(1232), more accurate data of this 
channel would give strong constraints on the range of cj and This nature of vr+p vr+vr+n 
can not be seen in other theoretical studies which have mainly focused on the vrA^ — > vrvrA^ reaction 
below = 400 MeV. 

The remarkable contribution of the process including the N^Tnr)^^.^^^^ vertex leads to the 
expectation that the ttN ttttN reaction becomes a source of information about the controversial 
scalar-isoscalar a meson. In the present work, the a meson is considered as a dynamical object 
generated by the rescattering mechanism of two-pion in the I = J = channel. We have obtained 
the following interesting results related to this meson: if we do not include the tttt rescattering effect 
in the N* N (inT)^^^^^ vertex, then (i) the large asymmetry in the vr'^vr'^ invariant mass distribution 
for ■K~p TT^TT^n can not be reproduced, and (ii) it is difficult to simultaneously describe the total 
cross section of all channels. The systematic analyses of all channels would provide an indication 
of the existence of a meson. 
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TABLE I: The matrix elements necessary to calculate the diagrams in Fig. |5] 





{N\jA\m 


{N\a\N) 


{N\n\N) 


{^\jv\m 










(A|j^|A) 







Here we mention the similarity about the role of A^*(1440) between the vrA'^ — > vrvrA'^ reac- 
tion and the A^A^ induced two-pion production reactions. The interference between A^*(1440) 
N(7nr )c~? and A^*(1440) Att is observed also in the pp ppir^ir^ and pn d{7nr)^^^ reac- 
tions UM ) actually plays a important role for explaining the data 



20|. In view of this 



similarity, we could also discuss the indications of a meson through those reactions. 

Finally, we comment on other baryon resonances related to the vrA^ ttttN reaction, which 
have been referred in several places. In the present work, we did not include the higher resonances 
such as A^*(1520) and A^*(1535). We have seen that, below the A^*(1440) mass-shell energy (i.e. 

< 480 MeV), the ttN ttttN data is almost saturated by the nucleon, A(1232) and iV*(1440). 
Therefore the higher resonances would become the relevant degrees of freedom at least above 

~ 480 MeV. Several discrepancies between the data and our results above T^^ = 500 MeV in the 
n^p vr^vr^p, n^p -K^ir^n and ii^p tt^-k^ti channels, would be cured by the consideration 
of such higher resonances. Anyway, we need further investigations on these resonances. 
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APPENDIX A: DETAILS OF THE MODEL IN NUCLEON-A SECTOR 

As mentioned in Subsec. Ill Al in this Appendix we explain some details of the model used in 
our previous study of vr ttttN Q] . To calculate the diagrams in Fig. [21 we need to evaluate 
seven matrix elements shown in Table d 
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1. Vector-isovector part 



The vector current matrix elements are 

{N{p')\j^^mNip)) = uip') 
for the nucleon, and 



-u{p) 



(Al) 



(A2) 



for the A^-A transition. In this appendix the isodoublet Dirac spinor for the nucleon is denoted as 

Based on the phenomenology of vector meson dominance (VMD) that at low energy the matrix 
elements of vector current jy are dominated by the p meson pole, we write the nucleon form factors 
as 

.2 



rin?p — t — impTpit) ' 



(A3) 



N 



(A4) 



where rup is the p meson mass and Ky is the isovector magnetic moment. The phenomenological 
width of the p meson is parameterized as 

3/2 



(A5) 



where T p is the total width at t = 



23] • As for the A- A transition form factors, we obtain 

2 



m?p-t- impTp{t) ' 



(A6) 



1 



F{^nt), 



(A7) 



rriM + rriA 

where we use the pA"A interaction HB4|) . The pNA coupling constant is denoted as fpNA, and 
_nds to the gauge coupling constant of the hidden local symmetry model tor the vector 
mesons [TJ, |75|. The other form factors Fy^{t) and Fy^{t) in Eq. HA2() are fixed to zero as long 
as we consider the Lagrangian ()B4|1 for the pNA interaction. 
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2. Axial-isovector part 



The matrix elements of are written as 



0,0 



{N{p'm^mN{p)) = u{p') [FUt)^, + Fl^{t)q,] j,-n{p), 
{A{p'm^mN{p)) = U''{p')[F^^Ht)9u, + F^^Ht)Q.7, 



+FAiii)QuQf. + iFlt{t)Q,<J^xQ 



NA, 



^ V2' 2 



u{p), 



(A(y)|jV(0)|A(p)) = U%p')[Fl,{t)g,y^,+Fl^{t)g,^q^ 

+Fii{t)qul^qx + FlS)q,q^qx] 75/"(i, W^iv)- 



(A8) 



(A9) 



(AlO) 



Same as in the case of N-N* and A-N* transitions caused by the axial current j^, some of the 
form factors in Eqs. HA8|) - HA10|) are exactly related to the renormalized coupling constants for the 
corresponding pion-baryon interaction, 



Jn 



(All) 



firNAit) 



[FZHt) + {m^ - m^)Fl^{t) + tFl^{t)\ , (A12) 



In 



2?7lA 



(A13) 



where we employ Eqs. (|Bip - HB3|) as the effective irNN, irNA and vrAA interactions. The other 
form factors (i.e. F^^{t), F^^{t), F^^{t) and F^^{t)) are fixed to zero. In our calculation all the 
form factors are taken as constants, and we eliminate F^^ and F^2 by using the PCAC hypothesis 
fnNA,7TAA{ml) ~ /ttATA.ttAaCO). 



3. Scalar-isoscalar part 

Since A(1232) dose not contribute to MsA, we only need the nucleon matrix element for the a 
current, 

{N{p')\amN{p)) = S{t)u{p')uip). (A14) 

According to the definition in Ref. Q], the form factor S{t) is equal to — cr,rAf(i)//7r'7i^, where 
(^nNit) is the pion- nucleon sigma term which becomes independent of t at tree level. 
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TABLE II: The value of constants used in the previous work. The mass and width of each particle or 
resonance are shown in MeV. 



Masses and widths 



(MeV) 



Parameters 





939 


U 


93 MeV 




138 


fp 


5.80" 


TOA 


1232 


KV 


3.71 


nip 


770 




45 MeV^ 


Ta 


120 




1.265 


Tp 


149 


pN 

pNA 

pNA 
^A,2 


5.67 X 10-3 MeV-i'= 

1.382'' 
-4.24 X 10-* MeV-i'^ 



"See p.33 iiv^ference Q 

''Reference 

''Reference 2£|. Note that the relation _F^2 = -2A^jv/mJ. 
"^Reference |n 

4. Pseudoscalar-isovector part 

The nucleoli matrix element of the vr current appearing in A^jr is written as 

{N{p')\t:%Q)\N(p)) = P{t)u{p')i^^T^u{p). 
The form factor P{t) are related to the renormalized ttNN coupling constant HAlOj) 

Pit) = ^ ( ^1 W(t), 



— t \ m-jr J 

where the pion pole contribution is taken into account. 

In Table. ^ we summarize the constants necessary to calculate the diagrams in Fig. |21 



(A15) 



(A16) 



APPENDIX B: PHENOMENOLOGY OF MESON-BARYON SYSTEM 

In this Appendix, we summarize the phenomenological Lagrangians used to estimate the form 
factors and mention a treatment of the finite width of baryon resonances. The Lagrangians of 
meson-baryon interaction are written as follows, 

fnNN 



L 



ttNN 



(Bl) 



C^NA = ^A-e,^(Zi)r (|, DNd^^TT- + H.c. , 



(B2) 
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C^AA = ^A^&a(i{Z2h^-f5p il l)e^5iZ2)A'dV, (B3) 

771 TT 



CpNA = i^A'^e.^(Z3)7.75/"(i, i)iV(9>^'^ - d^pn + H.c. , (B4) 



nip 



C^NN' = i^^N-f^^,T''N*dV + H.c, (B5) 



C^AN* = ^^^A^e,^(Z4)r (|, '^)N*dV + H.c, (B6) 

where and A^* are the isodoublet Dirac fields describing the nucleon and A^*(1440), respectively, 
and A'^ is the isoquadruplet Rarita-Schwinger field describing the A (1232). The second-rank 
Lorentz tensor 0^i, is defined by = (7^1, — ^(1 + 2Z)^^^y (where we take A = —1 

Q). The 

second term of this tensor vanishes if A(1232) is on the mass-shell (because of j^.U^lp) = 0), and 
so Z is called the off-shell parameter. In this paper we assume Zi = —1/2 (i = 1...4) for simplicity, 
i.e. 0^[y > dfiu- 

The A (1232) propagator is 



3(p^ — m^) 



(B7) 



rriA 

The nucleon and A^*(1440) propagators are expressed by the Dirac propagator which is familiar. 

In order to take account of the width of baryon resonances phenomenologically, we modify the 
denominator of the A(1232) and A^*(1440) propagators as — mj^ ^ — mj^ + imjiTfi^s), where 
i? = A or N* . The widths rA(s) and Tn*{s) are taken as 



/s |q(mA)|^ 



TAis) = rA-^ :^Z ' (^{V~s -TUN- m^), (B8) 



and 



rN'{s)=TN'-l^^^^e{V^-mN-m^), (B9) 



respectively. Here q = q(-v/s) is the pion spatial momentum in the center of mass ttN system with 
the total energy ^/s. In this paper we use m^* = 1440 MeV and Tat* = 350 MeV (as for the value 
of mA and Fa, see Table. HH). 
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Finally we list the formulae for the isospin matrices used in our calculation (see e.g. appendix 
AinRef. fl), 

P\1\)I\1\) = 5'^'-\t'^t\ (BIO) 

c; 1 9 1 

J- (2' 2)-' (2 ' 2)^ (2' 2) = ~ 6 3 6 ^ 
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